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Abstract 

Here diffusion and Ware pinch are analyzed as opposed effects for plasma confine- 
ment, when instabilities are not considered. In this work it is studied the equilibrium 
inductive electric field where both effects annul each other in the sense that the av- 
erage normal velocity is zero, that is, marginal velocity confinement is reached. The 
critical electric field defined in that way is studied for different values of elliptic 
elongation, Shafranov shift and triangularity. A short complementary analysis is 
also performed of the variation of the poloidal magnetic field along a magnetic line. 
Magnetohydrodynamic transport theory in the collisional regime is used as in recent 
publications. Axisymmetric and up-down symmetry are assumed. 



I Introduction 

The H-mode is characterized by the suppression of anomalous transport in tokamaks, 
because of low plasma turbulence induces by internal barriers [1-3]. As result neoclassical 
transport calculations becomes very important in this mode. Diffusion in the collisional 
regime depends of the pressure gradient at the 95% surface adjacent to the scrape of 
layer or SOL, as well as the inductive electric field E^^. The diffusion due to the gradient 
pressure is opposed to Ware the pinch effect due to E^. 

In previous papers was shown that neoclassical diffusion can be treated with great 
simplicity in the case of arbitrary plasma configuration, using a new kind of tokamaks 
coordinates described there [4-6]. However, suitable direct numerical calculations for dif- 
ferent values of elongation, Shafranov shift and triangularity have not been presented 
until now. Some previous results on this theme using these coordinates are incomplete, 
and they are not using the right parameters for the tokamaks in operation at present. 
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Here calculations are presented in a different way, since we look for the values of E^p in 
the marginal velocity confinement, that is, when the average velocity on 95 % surface is 
zero the results corresponding to marginal confinement flux, that is, the transition from 
the outgoing to ingoing flux, could be more interested, but much more difficult to calculate 
and a suitable and simple treatment describing this process for any plasma configuration 
seem that they have not been publishing until now. 

In the calculations now presented there are first a suitable normalization procedure 
absent in previous calculations as well as an adequate selection of tokamak parameters. 
The normalization used here allows us to get results, which can be useful for a diversity 
of different tokamak plasma configurations, with different values of ellipticity, Shafranov 
shift and triangularity. 



II Theoretical Treatment 

The collisonal transport treatment presented in previous paper for toroidal axisymmetric 
plasmas can be written in a more convenient way using dimensionless integrals and vari- 
ables, as 

^ = ^ (02^ = T^TWTR^ + ^VlR /ff ^2 (^3 - y )] p(/7 - ^ + 71/5) , (1) 

where v is the dimensionless normal velocity derived from the velocity < v > along the 
magnetic surface normal, and E^i is a dimensionless electric field defined as 

E^ 

^<P = r?x (^p^ ■ (2) 

Her B^pi, Bpi, E^i and (^)i are respectively the major radius, toroidal and poloidal 
magnetic field, inductive electric field and pressure gradient at the mid-plane external 
point Aq of the magnetic cross section. The plasma resistivity are r}± and r/y in the di- 
rections perpendicular and parallel to the magnetic field lines and Rc is the minor axis 
radius. The dimensionless quantity 71 is the ratio 71 = Bpi/B^i between poloidal and 
toroidal magnetic field at point Aq. On the other hand the new dimensionless integrals 
li, i — 1 to 7, are defined as 

f R(s) ds 

' f Rl ' ^ ^ 
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R^(s) ii(s) ds 

Rt (1 + 7? KsY) 
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r I /^(g) Rjs) ds 

5 - f ^ir; — : — 2 — ' 



ds 



fi{s)R{s) ' 



(9) 



/ ^ , (10) 

where all the integrals are around a magnetic surface and //(s) is a function, depending 
of the curvature Ka- of the orthogonal line family, giving by 

= exp { — K-a ds) . (11) 

This results are obtained using MHD equations and assuming toroidal axisymmetry. 

In order to get numerical as well as analytic results it is useful to express the family 
of magnetic cross sections by the equations 

R/Rc = 1 + X[(E - l)cose + Tcos{2e) - A] , (12) 

z/R^ = \[{E - l)sin9 + Tsin{2e)] , (13) 

where E, T and A are respectively ellipticity, triangularity and Shafranov shift distortions. 
The parameter A in this equations labels each magnetic surface. The previous equations 
are a generalization of the equations presented by Roach, et al[7] . The quantities E, T 
and A are dimensionless, however for the analysis and calculations are more useful the 
new dimensionless quantities E, T, A, R and z, defined respectively as 

~ A ~ T ~ , ^ ~ R z 

^ = ^ ^ r = — ; A = A(£-1) ; ; i = - . (14) 
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The well know Shafranov shift Ashaf is connected to the previous parameters by 

Ashaf = {A-f)a , (15) 

where 2a is the size of the 95 % magnetic surface measurement at the mid-plane, or in 
different words, a is the horizontal half-width of the plasma. The equations of the cross 
section magnetic lines will be now 

^ = 1 + ~X[cose + fcos{2e) - A] , (16) 
^ + 1 



z 



—— X[sine + Tsin{2e)] . (17) 
E — 1 



Denoting by the value of the parameter A generating the 95 % magnetic surface with 
general coordinates Za, then 

Ra ^ 1 + Xa[cose + fcos{29) - A] , (18) 

^a= {^^) k[sine + fsin{2e)] . (19) 
The largest and smallest values of Ra are respectively Rai and Ra2, and its values are 

^al=^a(^ = 0) = l + A„[l+r-A] , (20) 
^a2=^a(^ = 7r) = l + Aj-l+T'-A] . (21) 

The radius Rao of the center of the plasma and plasma size a will be respectively 

Rao^ ^"' t ^°° -l + A„[l + f-A] , (22) 



a _ Ral — RaO _ T /r,g\ 

R;- 2 ■ ^^^^ 

It seems also convenient to connect that parameters with the aspect ratio A, such that 

, Rao (1 + A„(f-A) 1 . T ^ 

A = — = ^ ^ ^ ; ^ A + A-T . (24) 

« Xa Xa 

As in our previous papers, b is the maximum value of Za- The elliptic elongation K 
and triangularity are connected to our previous parameters as 



E + 1 /3 1 
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where g and h are 

5 = 32r'(§i|/ . A=i(yrTi-l) (26) 

and 

where Ram is the radius of the point with maximum z. 



(27) 



III Results 

The dimensionless variables defined previously simplifies the computation, because quan- 
tities we need for the calculation are: the ratio between poloidal and toroidal magnetic 
fields 7i, the horizontal half- width of the plasma a, the aspect ratio A, the ratio r)\\/r)±, 
the Shafranov shift, the ellipticity K and triangularity S. Prom K and S, the values of 
E and T are determined. The Shafranov shift, the value of a and the calculated value T 
together allow the calculation of A and Aq. 

In this way the family of magnetic surfaces can be determined. Following the proce- 
dures described in previous works the family of orthogonal lines can be also determined, 
which allows to obtain the curvature function and the function All the integrals 
needed to obtain the normal dimensionless velocity v can also be performed once the value 
of 7i is given. The velocity v can also be obtained if the ratio f]\\/r]± is given, and it is 
a linear function of the dimensionless toroidal electric field E^i. The intersection of that 
value with the axis of abscissa allows to determine the critical dimensionless electric field 
Eipicrit, for marginal velocity confinement. This field will be later determine for different 
values of ellipticity K, dimensionless Shafranov shift i^shafrancm/ and triangularity 5. 

Following the above described procedure, the ellipticity K and triangularity 5, in Fig- 
ure 1, are given as i^' = 1.76 and 5 = 0.25. These correspond to values of E and T: 
E = i, T = 0.3. The real procedure we use was a little different. We first select values 
of E and T, in such way K and S become about their values in JET tokamak[8]. This 
procedure is simpler for us, and it is the same idea. 

In Figure 1, several cross sections magnetic lines have been drawn for different values 
of A in the interval from zero to A^, giving in Eq.(24). The characteristic values for that 
figure are: elipticity K = 1.76; relative Shafranov shift Ashaf = ^shaf/a = 0.3; triangu- 
larity S — 0.25; horizontal half width a — 1.12 m ( 95 % surface ); aspect ratio A — 2.5 
and minor magnetic axis ratio Rc — S m. From the previous data, the following param- 
eters are determined: elliptic dispersion E = A, sec Eq.(25)to (27), relative triangularity 
dispersion T — T/{E — 1) = 0.3; relative Shafranov dispersion A = ^shaf + T = 0.9; 
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relative A-parameter at the 95 % surface \a = {A + A — T) — 0.1075; the outward 
radius Ra^, Rl = R{0 = 0) = Rao = -Rc[l + Ki^ + f - A)] = 3.38 m; the inward 
radius R2 = R{6 = vr) = Rc{—^ + A(i(l + T — A)] = 1.45 m; center plasma radius 
Rq = {R1+R2) /I = 2.42 m and radius at the maximum z, {{dz/dR)^^ = 0),Rm = 2.18 m. 

In Figure 2, the function /^(s) is show as a function of 9. This function allows to deter- 
mine poloidal field around a magnetic surface, once the value Bpi at the outward point Aq 
is measured or determined. In order to find ii{9) , it is necessary to determine the curvature 
of the family of orthogonal lines. The procedure has been explained elsewhere[4]. The 
function fi{6) appears in most of the integrals needed to determine the average normal 
velocity < v > first calculated by Pfirsch-Schliiter for cross-sections of circular magnetic 
surfaces[9]. The function to be used in this work is the central function, plain line. How- 
ever, two other //-functions are also shown to illustrate that function. The upper curve 
correspond to the case where the triangularity and Shafranov shift are zero. The lower 
curve is just the case of zero triangularity, but the same Shafranov shift than in the main 
curve ( central one), where the triangularity is 5 = 0.25, as in Figure 1. 

Since the function /x(s) shown essentially the behavior of RBp, the upper curve illus- 
trate that this product is constant at the inward and outward points, when there is not 
triangularity and Shafranov shift. However, the value of RBp, decreases at the uppest 
and lowest points because of the ellipticity elongation K = 1.76. The Shafranov shift 
modifies this pattern in the way that the values at the inward area are almost constant, 
and very low compared with those at the outward point Aq. Introducing triangularity as 
in the case of the central line does not modify the pattern, but the values at the flat part 
of the curve are not so low. 

In Figure 3, the dimensionlcss average normal velocity is shown as a function of the 
dimensionless toroidal electric field E^i at point ^40, the line here found is the straight 
line, because of the normalization here used using the value of {r]±/B^){dp/da)i at the 
point Aq. The intersection of that line with the abscissa define the critical dimensionless 
toroidal field E^i for marginal velocity confinement. This critical electric field is show as 
a function of the elliptic elongation K, where all the other variables are kept fixed at the 
values shown in Figure 1, that is, 6 = 0.25 and A = 0.9. To determine the integrals in or- 
der to draw Figure 3 and 4, the value of 71 is chosen as 0.3, which correspond for instance 
to toroidal magnetic field B^i — 5 Teslas and poloidal magnetic field Bpi — 1.5 Teslas. 
The ratio of perpendicular to parallel resistivity, r/y/r^x has been considered as 1.97 as in 
page 669 of Ref.(8). As a way of completion Tg is taken as 5 keV. 

The Figure 3 shown that the outward velocity due to diffusion is opposite by the inward 
effect due to the inductive electric fields E^. After a critical value E^,, the characteristic 
Ware pinch effect [10] becomes more important, and an average inward velocity < v > is 
produced, in such a way that the plasma appears confined as long as instabilities arc not 
considered. The critical dimensionless electric field increase with the elliptic elongation K 
if all the parameters are kept fixed, however, the variation is not so large, from 0.8 a factor 
to about 1.1 as it is shown in Figure 4. In this figure has been normalized with a second 



6 



procedure. First a normalized critical toroidal electric field is selected as reference and 
denoted by E^^i cHt. ref.i which in this case corresponds to that given in Fig. 3, where the 
characteristic values of the parameters above given. The value of this critical dimension- 
less toroidal electric field is 29.8106, and correspond also to the horizontal line through 
one, which it is also show by a dot line. This kind of normalization is also performed in 
Figure 5 and 6. The toroidal electric fields with the second normahzation explain above 
are denoted by E^p\ crit.- 

The changes due to Shafranov shift with all the other parameters fixed, are more sig- 
nificative as illustrates the Figure 5, where the critical electric field could be one forth 
or 3 times the value of that shown in Figure 3. Furthermore the curve in Figure 4 is 
almost linear, but not that in Figure 5, which seems somewhat as a parabola. Finally, in 
Figure 6, the change of E^crit with triangularity are shown. There also crit. changes 
strongly with S. More important than this it is that the changes are very significative for 
low values of S, and E^i crit. could be lower by a factor 5, and with values of triangularity 
no so large, as 6 = 0.25. Here also the curve is strongly not linear, however the concavity 
of the curve is opposite to that in Figure 6. 



IV Conclusion 

In most tokamak operation the inductive magnetic field effect exceeds the plasma diffusion 
effect and the Ware pinch effect contract the toroidal plasma column. Here the critical 
point where both effect becomes almost equals is studied. This equilibrium situation is 
consider as that where the average normal velocity becomes null or void, which will be 
defined as the marginal velocity confinement. A suitable normalization procedure allows 
to extend our analysis to a large amount of different situations in tokamak plasma config- 
urations. The critical toroidal electric field changes very little with the ellipticity of the 
plasma. However, the changes are very strong with the Shafranov shift and triangularity. 
The curves in those last cases seems parabolas, but with opposite curvature. Very large 
changes for small values of triangularity have been found, producing changes with a factor 
5 for modest triangularity values as S — 0.25. 
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Figure 1: Magnetic flux surfaces showing the characteristic radius, triangularity and main 
point, with the parameters values given in the text. 
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Figure 2: Characteristic exponential factor around a magnetic surface for the main val- 
ues of the reference surfaces, plain line, and two other situations. No triangularity but 
Shafranov shift, dashed line, and neither Shafranov and triangularity, point line. The 
parameters values of the reference curve are: E — 4, A — 0.9, A — 2.5, T — 0.3, 71 = 0.3 
and Rc = 3 m 
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Figure 3: Normalized velocity versus normalized toroidal electric field in the outward 
point for the parameters values described as reference parameters in the text and for the 
outward. 
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Figure 4: Twofold normalized critical electric field E^i cru. as a function of elliptic elon- 
gation with the triangularity and Shafranov shift equal to those in the reference curves, 
Fig. 1 and Fig. 3 
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Figure 5: Twofold normalized critical electric field E^^i crit. as a function of dimensionless 
Shafranov shift with elliptic elongation and Shafranov shift given by the reference curve 



13 



9- 



2 - 




0.1 



0.2 



0.3 



0.4 



0.5 



0.6 



Figure 6: Twofold normalized critical electric field E^i crit. as a function of triangularity 
with ellipticity elongation and Shafranov shift given by the reference curve. 
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